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Abstract 

We estimate the most important corrections to the branching ratios for the 
invisible decays of quarkonium states, arising from possible extensions of the 
Standard Model. Among the possibilities considered are the presence of extra 
Z-bosons, minimal supersymmetric extensions of the Standard Model with 
R-parity violation and decays into Goldstinos. Prospects of detecting these 
corrections at existing and future B-factories and r-charm factories are dis- 
cussed. 
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B-meson factories under construction at KEK and at SLAC [1] can deliver 10 8 T's and 
T"s per year, when they are tuned to run at resonance. Coupled with what is available at 
CLEO III, these facilities will make it possible for the first time to study in detail invisible 
decays of these resonances. Such decay modes can be studied by observing the decay T' — > 
T + 2n's at resonance and tagging on the invariant mass of the dipion system at the T-mass. 

Invisible decays of heavy quarkonium states offer a window into what may lie beyond the 
Standard Model. The reason is that apart from neutrinos, the minimal Standard Model 
predicts no other channels that these states can decay into. The associated rates for these 
decays can be computed precisely, and so any observed departure can furnish hints of struc- 
tures over and above those in the Standard Model. As we will show below, a more complete 
test of what these structures might be would require similar measurements on the J/ty- 
invisible decay widths, something which could be achieved at the r-charm factory in the 
future. 

In what follows, we concentrate on what we believe to be the largest among these effects. 
In particular, we will estimate the branching ratios into neutrinos in the presence of extra 
Z-bosons, and R-parity violating effects in supersymmetrized Standard Models, and finally 
decays into Goldstinos. 

To begin, we present the Standard Model prediction for the branching ratio of the in- 
visible decays of T and J/ty and their observed decays into electron-positron pairs. To 
our knowledge the detail formulae have not been given before. Within the Standard Model 
the invisible mode consists solely of decays into three types of neutrino-antineutrino pairs. 
Neglecting polarization effects and taking into account e+e - production through a photon 
only we get 
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with G and a being the Fermi and the fine structure constants respectively. Mx^m are 
masses for the T and J/\& states. These formulas are expected to be correct up to about 
2-3%. The major sources of theoretical uncertainty involved in (1) and (2) can be listed as 
follows: 



1. corrections to the T and J/ty wave functions, including QCD corrections and polar- 
ization effects, 

2. e+e - production through Z, 

3. electroweak radiative corrections, 

4. corrections due to the Higgs boson. 

The first correction is not expected to be significant since it similarly modifies the neutrino 
and electron decay widths, and so the leading order cancels in the branching ratio. The 
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second correction will introduce a factor O.lMg/Mf , and so is not expected to be significant 
either. 

The third correction can be substantial, and thus needs to be considered in a little more 
detail. In the numerical calculation quoted above, we have included the MS'-running of the 
coupling constants in a, G, as well as in sin 2 9. We have chosen the Z-mass, G, and a(M z ) 
as the input parameters and run the other parameters from their values at the Z-mass 
down to the low energies where the quarkonia reside. These are by far the most significant 
contributions. To the accuracy we are working in, of the order of 2-3%, we may ignore 
threshold effects. Box-type diagrams with double Z (W + ,W~) or photon emission which 
may be expected to modify the tree level results turn out to be negligible. The largest of 
these with two internal photon lines vanish due to charge- conjugation symmetry since both 
T and J/\& are C-odd. The contribution of the remaining box diagram with virtual Z and 
7 also makes only a negligible correction to the tree level Z-mediated e + e~-production and 
can be discarded as well. Finally, the Higgs correction is vanishingly small simply because 
we are dealing with virtually massless particles in the final state. 

As a result, branching ratios of (1) and (2) are theoretically clean and thus offer a rare 
opportunity to search for physics beyond the Standard Model if they are relatively large. 
Among the most viable candidates are supersymmetry (SUSY) and Grand Unified Theories 
(GUTs). In the case of spontaneously broken SUSY the decays of T and J/ty will produce 
invisible Goldstinos (or light gravitinos) in addition to the neutrino models. All other par- 
ticles in the final state are prohibited by kinematical considerations (see [2]). Therefore, 
measuring the invisible width and comparing it with the SM prediction could provide in- 
formation about new physics. In the case of SUSY models, we can gain information on 
the SUSY-breaking scale as well as about the existence of R-parity violating terms in the 
superpotential. 

We first consider the case of decay into light gravitinos denoted by g. The general structure 
of the Goldstino interactions is known, and given for example in [3]. 

T and J/\& can decay into Goldstinos via virtual Z, 7 in the s-channel and via exchanges 
of b, c-squarks in the t-channel. Neglecting the Goldstino exchange effects, the corresponding 
rates for the J/^f are calculated to be 
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where v 2 = v \ + v \ ~ (174 GeV) 2 , v 12 = ($J 2 ), tan/5 = v 2 /vi, m c (rh c ) is the mass of the 
c-quark (squark) and F is the SUSY-breaking scale. The photon channel is suppressed as 
compared to the Z and c ones since the photon-goldstino coupling containes higher powers of 
F and can be neglected in the leading order approximation. Unfortunately, for a reasonable 
choice for the value of F ~ 1 TeV[3], these rates are extremely small and far beyond the 
experimental capabilities. Doing the same calculation for the T does not improve matters 
for gravitino decay modes. The above considerations rule out light gravitinos as candidates 
for the invisible decays. 

Nevertheless, supersymmetry still can affect the invisible width because of R-parity break- 
ing processes. Such processes affect neutrino decay modes through squark exchange. In 
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terms of superflelds, the relevant interactions are generated by a A^L^Qi-D^ term in the 
superpotential. Here i,j, k are the generation indices and we have suppressed SU(2) indices 
[4]. Expressed in terms of component fields the interaction Lagrangian takes the form 

— C R = Xij k [v iL d kR dj L + dj L d kR v iL + d* kR v c iL dj L 

- (e iL d kR u jL + u jL d kR e iL + dl R (^ iL )u jL )\ + H.c. (5) 

where we have neglected mixing among generations. From Eq. (6) it can be seen that the 
terms leading to neutrino final states involve the down quarks only, and the J/\l/ width will 
not be affected. It is noteworthy that one needs to take into account non-SM corrections 
to the neutrino widths only, since their contributions to e + e~ production are far less than 
0.1% and lead to higher order corrections to the branching ratio. Neglecting possible squark 
mixings and quark mixings we get 
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One notices that R-breaking contributions add coherently to the SM result and reduce the 
width. The present experimental constraints on A- 33 are very loose [5] and SUSY corrections 
may turn out to be quite significant for this reaction. For a SUSY mass of 100 GeV, 
A' 133 < 0.002 from v e mass calculation, A 233 < 0.4 and A 333 < 0.26 from the ratio of hadronic 
to leptonic widths at the Z-pole [6]. We note that Eq. (6) also includes an incoherent piece 
coming from decays into neutrinos of different flavors. 

We display in Fig. 1 the branching ratio of Eq. (1) as a function of the parameter 

3 / 2 

x(SUSY) = ^2 J^<3|L where we have set m~ b = m~ b = m~ b = 100GeV/c 2 for definiteness. 

i=i m ' b " L 

It is seen that corrections to the SM result as large as 30% are possible for a range of Im- 
parity violating couplings. We have cut off the x values at the maximally allowed value of 
z max = 2.2 x 10 -5 GeV~ 2 as dictated by experimental bounds only the A"s. Similarly the 
sensitivity of the invisible T decay to the Z' mass is given in Fig. lb where now X = M\jM\,. 

We next consider another unconventional contribution to T {T,J/^/ — > vv) - an extra 
neutral gauge boson, Z\ which provides an additional annihilation channel. Z' bosons 
appear as remnants of a higher symmetry at large energies and have to be sufficiently 
massive in order to fit current experimental limits. We will concentrate mainly on the 
superstring-inspired Eq grand unification model [10,7] and left-right symmetric models [8]. 
In the models under consideration, the Z' will correspond to an extra U(l) (for Eq) or to a 
neutral component of SU(2) R (for left-right models). 

The phenomenology of an extra Z 1 is highly model dependent [9]. For the case of Eq 
this depends on the breaking of Eq to the SM group. The details are beyond the scope of 
this paper. We are only interested in probing Z 1 with Eq. (1) and (2). To this end, we 
first note that we can neglect the Z — Z' mixing because precision measurements at LEP2 
put a bound of < 0.0025 for such mixings [9]. The fermion couplings and hypercharges are 
uniquely determined by the way the Standard Model is embedded in Eq [10]. The neutral 
current process mediated by the Z' boson involving the fermion / is given by 
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Cnc = g E {Y' fL f L Yf L + Y' fR f R ^f R )Z^ . (7) 

where qe = ^J^g tan. 9 W , g is the SU(2) gauge coupling, and Yj-'s are the hypercharge of 
the fermions of a given chirality. The relevant charges for the superstring-inspired Eq are: 
Y' , =\J\, Y' =-\Jl, Y' = -\J\ and Yl = k/f. The result is 

c l,°l 3 V 5 ' cr 3 V 5 ' v 6 V 5 or 6 V 5 

T ( 2fi 2 + 
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where ft = -1 + f sin 2 9 W + 4cos2 J™ M % 9 \y' v (JI l + F 6 ' ). Interestingly, the SM result for 

J/\I/ width is unaffected by the extra Eq-Z' boson since Y' + Y^' = 0. In contrast with R- 
breaking processes, the Z' increases the T width. Hence, if both of these new physics sources 
are present, destructive interference between these non-SM corrections can take place if the 
parameters are favorable. This would be fortuitous though not impossible. 

For the case of the left-right model, both (1) and (2) undergo a certain modification. 
Assuming -^f- <C 1, the mixing angle between Z and Z' can be expressed as [8] 

M\ 
Ml 

Then 



~ v/cos2^ — f . (9) 
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One readily obtains the corresponding expression for the T decay width from (8) with 

n => (-1 + ^sin 2 9 w )A . (11) 

In this case the widths for both quarkonia will decrease as compared to their SM values. 

In Fig. 2 we show the sensitivity of the T decay to the extra Z boson in the left-right 
symmetric model with x(LR) = M|/M§,. Very similar behaviour for the branching ratio of 
J/ip decay is displayed in Fig. 2b. Evidently, the latter resonance will not be suitable for 
probing this class of models. 

We have argued that because of the theorectically clean nature of the decays a careful 
measurement of the invisible widths of the heavy quarkonium states can therefore yield 
constraints on a variety of physics beyond the Standard Model. It is especially useful for 
studying R-parity breaking terms of the third generation in SUSY models. It is also sensitive 
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to extra Z bosons of GUTs or the left-right symmetric model. We also showed that the 
invisible channel will have to be due to the light neutrinos since the other possibility of light 
gravitinos will not be significant. Neutralinos will also not contribute since a lower bound 
of 40 GeV has already been established by LEP measurements [11]. If a deviation from the 
SM value is found in the T decay a similar measurement for the J/ip can shed light on the 
source of new physics. 

The work is supported in part by the US Department of Energy under Grant No. DE- 
FG05-92ER40709-A005, and the Natural Science and Research Council of Canada. 
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I. FIGURE CAPTIONS 



Fig. la The branching ratio of T — > uu/T — > e + e~ as a function of the parameter 
x(SUSY) = Y,i I -M 33 I 2 / m i m R-parity violating MSSM. A' and the m~ b are given in 
the text. 

Fig. lb The branching ratio T — > i/V/T — > e + e~ as a function of rr = M§/Mf, for the extra 
Z'-boson in E(6) models. 

Fig. 2a The branching ratio T — > uV/T — > e + e _ as a function of x{LR) = M^/M' z 2 in the 
left-right symmetric model. 

Fig. 2b The branching ratio of J/^ — > i/X// J '/\& — > e + e~ as a function of x(LR). 
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